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Abstract. Noise, through its interaction with the nonlinearity of the living sys- 
tems, can give rise to counter-intuitive phenomena such as stochastic resonance, 
noise-delayed extinction, temporal oscillations, and spatial patterns. In this paper 
we briefly review the noise-induced effects in three different ecosystems; (i) two 
competing species; (ii) three interacting species, one predator and two preys, and 
(iii) N-interacting species. The transient dynamics of these ecosystems are ana- 
lyzed through generalized Lotka-Volterra equations in the presence of multiplicative 
noise, which models the interaction between the species and the environment. The 
interaction parameter between the species is random in cases (i) and (iii), and a 
periodical function, which accounts for the environmental temperature, in case (ii). 
We find noise-induced phenomena such as quasi-deterministic oscillations, stochastic 
resonance, noise-delayed extinction, and noise-induced pattern formation with non- 
monotonic behaviors of patterns areas and of the density correlation as a function 
of the multiplicative noise intensity. The asymptotic behavior of the time average of 
the j"" population when the ecosystem is composed of a great number of interacting 
species is obtained and the effect of the noise on the asymptotic probability distri- 
butions of the populations is discussed. 



1. Introduction. In recent years several theoretical investigations have been done 
on noise-induced effects in population dynamics CD|-IH|- In particular, the problem of 
the stability of complex ecological systems in the presence of noise has been widely 
discussed P|. New counterintuitive phenomena, such as stochastic resonance ^1 
[TT] . noise enhanced stability ^21 and noise delayed extinction 0^1^], can appear 
because of the presence of noise in living systems, whose dynamics is nonlinear. 
The interaction between noise and nonlinear determinism in ecological dynamics 
adds an extra level of complexity compared with the largely stochastic dynamics 
of, say, economic systems or the largely deterministic dynamics of many physical 
and chemical processes QHI- Ecological systems are open systems in which the 
interaction between the component parts is nonlinear and the interaction with 
the environment is noisy. This intrinsic nonlinearity can give rise to the complex 
behavior of the system, which becomes very sensitive to initial conditions, various 
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deterministic external perturbations, and to fluctuations always present in nature. 
The comprehension of noise's role in the dynamics of nonlinear systems plays a key 
aspect in the efforts devoted to understand and then to model so-called complex 
ecosystems. One approach to understanding the complexity is to start with a 
conceptually simple view of the system in order to catch the phenomena of interest 
and then to add details that introduce new levels of complexity E] ■ In general 
the effects of small perturbations and noise, which are ubiquitous in real systems, 
can be quite difhcult to predict and often yield counterintuitive behavior. Even 
low-dimensional systems exhibit a huge variety of noise-driven phenomena, ranging 
from a less ordered to a more ordered system dynamics. 

In the past, the study of deterministic mathematical models of ecosystems has 
clearly revealed a large variety of phenomena, ranging from deterministic chaos to 
the presence of a spatial organization. These models, however, do not account for 
the effects of noise despite the facts that it is always present in actual population 
dynamics and that it arises from different sources, such as the intrinsic stochasticity 
associated with the random variability of the environment. Frequently, its effects 
have been assumed to be only a source of disorder. Recently researchers have shown 
a growing interest in a deeper understanding of the eflfects of fluctuations in bio- 
logical systems ranging from neuroscience to biological evolution and to population 
dynamics [T]-|H1, [I1]-EI]. 

In addition, analyses of experimental data of population dynamics frequently 
need to consider spatial heterogeneity. Characterizing the resultant spatio-temporal 
patterns is, perhaps, the major challenge for ecological time series analysis and for 
dynamics modeling. To describe complex ecosystems, it is therefore fundamental 
to understand the interplay between noise, periodic and random modulations of 
some environment parameters, and the intrinsic nonlincarity of simple models of 
ecosystems and to understand spatio-temporal dynamics P [7j.[^-P()]. 

The principal aim of this work is to review some recent results obtained for 
systems described in term of a generalized Lotka-Volterra model including a term 
of multiplicative noise A constructive role of the noise is observed. It con- 
tributes to producing: (a) quasi-periodic oscillations and stochastic resonance in 
the presence of a driving force; (b) noise-delayed extinction, (i.e. a nonmonotonic 
behavior of the average extinction time of one of the two species as a function of 
the noise intensity); and (c) nonmonotonic behavior of the pattern formation, the 
density correlation and pattern areas as a function of noise intensity. We analyze 
three different ecosystems described by the formalism of the Lotka-Volterra equa- 
tions. The first ecosystem is comprised of two competing species in the presence 
of two noise sources: a multiplicative noise which affects directly the dynamics of 
the species, and an additive noise responsible for the random behavior of the in- 
teraction parameter between the species. We obtain quasi-periodic oscillations of 
two species densities, stochastic resonance (SR) and noise-delayed extinction. We 
also investigate the system using multiplicative colored noise with different values 
of the correlation time Tc- The effect of the correlated noise is to shift the peak of 
the signal-to-noise ratio (SNR), which is the signature of the SR phenomenon. For 
this ecosystem we also analyzed the spatial effects by considering a discrete time 
evolution model of the Lotka-Volterra equations with diffusive terms, namely a cou- 
pled map lattice (CML), and we analyzed the spatio-temporal patterns of the two 
species induced by the noise. In the second ecosystem, comprised of three interact- 
ing species, namely one predator and two preys, we analyzed the spatio-temporal 
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behavior of the species densities. Wc find: (a) noise-induced pattern formation 
in the coexistence regime, which depends on the initial conditions, (b) oscillating 
behavior of the site correlation coefficient with an alternation between coexistence 
and exclusion regime, (c) nonmonotonic behavior of the pattern area as a func- 
tion of noise intensity. Finally we consider a system comprised of many interacting 
species. The analytical resolution of the Lotka-Volterra equations is more difficult 
in the presence of a large number of species. Nevertheless some analytical approx- 
imations for the mean field interaction between the species as well as numerical 
simulations give some insight into the behavior of complex ecosystems |H1 EH . 
For a large number of interacting species, it is reasonable, as a phenomenological 
approach, to choose the growth parameter and the interaction parameter at ran- 
dom from given probability distributions. Within this type of representation, the 
dynamics of coevolving species can be characterized by statistical properties over 
different realizations of parameter sets. Though the generalized Lotka-Volterra 
model has been explored in detail |^, it seems that a full characterization, either 
deterministic or statistical, of the conditions under which a population becomes 
extinguished or survives in the competition process, has not been achieved fZP] I30j . 
In this last ecosystem, two types of interaction between the species have been con- 
sidered: (a) mean field interaction, and (b) random interaction. We focused on the 
statistical properties of the i*'' population, obtaining the asymptotic behavior of 
the time integral and the distributions both of the population and the local field, 
which is the interaction of all species on the population. By introducing an 
approximation for the time integral of the average species concentration M{t) we 
obtained analytical results for the transient behavior and the asymptotic statistical 
properties of the time average of the i*^ population. 

2. Two competing species. Time evolution of two competing species is obtained 
within the formalism of the Lotka-Volterra equations |31| in the presence of a mul- 
tiplicative noise 



— = fJ.ix{ai - X - (3i{t)y) + xCj;{t) (1) 
at 

dy 

— = tJ-2yia2~y- P2it)x) +y^y{t), (2) 

where (t) and (t) are statistically independent Gaussian white noises with zero 
mean and correlation function {^i{t)$,j(t')) = ad{t — t')5ij = x,y). It is known 
that the real biological systems are affected by random interactions, because of the 
presence of environmental fluctuations. The noise and some other deterministic 
periodical driving force present in the ecosystems, such as the temperature, con- 
tribute to determine also the dynamics of (3, the interaction parameter between the 
species. For /3 < 1 a coexistence regime takes place (that is, both species survives), 
while for /? > 1 an exclusion regime is established (that is, one of the two species 
vanishes after a certain time). Coexistence and exclusion of one of the two species 
correspond to stable states of the Lotka-Volterra's deterministic model The 
change in the competition rate between exclusion and coexistence occurs randomly 
because of the coupling between the limiting resources and the noisy environment. 
A random variation of limiting resources produces a random competition between 
the species. The noise therefore, together with the periodic force, determines the 
crossing from a dynamical regime (/3 < 1, coexistence) to the other one (/3 > 1, 
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Figure 1. The bistable potential U{/3) of the interaction parame- 
ter (3{t). The potential U{(3) is centered on (3 = 0.99. The param- 
eters of the potential are h = 6.25 • 10"'^, rj = 0.05, p = —0.01. 

exclusion) [3 . To describe this continuous and noisy behavior of the interac- 
tion parameter P(t) we consider a stochastic differential equation with a bistable 
potential and a periodical driving force 

~dr ^ df + Tcos (wot ) + C/3 ) , (3) 

where U{13) is a bistable potential (see Fig. 1) 

[/(/?) = h{(3 - (1 + p))V^^ - 2h{(3 - (1 + p) f/rj', (4) 
and h is the height of the potential barrier. The periodic term takes into account 
for the environment temperature variation. Here 7 = 10^^ and loq/ (27r) = 10"'^. In 
Equation |j2J) ^i3{t) is a Gaussian white noise with the usual statistical properties: 
(C/3(t)) = E^iid {^0{t)^fi{t')) = ap6{t — t'). Due to the shape of U{(3) it is reasonable 
to expect a coexistence regime for /3(0) < 1, when deterministic case {£,p{t) = 0) is 
considered. 

2.1. Stochastic resonance. First, we investigate the effect of the noise on the 
time behavior of the species. Since the dynamics of the species strongly depends 
on the value of the interaction parameter, we initially analyze the time evolution of 
(3{t) for different levels of the additive noise ap. To obtain the time series for the two 
species, we set in Equation (jJl ai = 0:2 = a, (3i[t) = (32{t) = I3{t). Depending on 
the value of the multiplicative noise intensity we obtain: (i) a periodical behavior 
of /3(t) in the coexistence region (see Fig. 2a); (ii) the same behavior of Fig. 
2a, slightly perturbed by the noise (see Fig. 2b); (iii) a quasi-periodical behavior 
of the interaction parameter jumping between the two values (3 = 0.94 < 1 and 
/3 = 1.04 > 1, respectively corresponding to left side well (coexistence regime) and 
right side well (exclusion regime) of the potential shown in Fig. 1; and finally 
(iv) a loss of coherence and a dynamical behavior strongly controlled by the noise 
(Fig. 2d). We note in Fig. 2c synchronization of noise with driving periodical 
force ^1 m , the typical signature of stochastic resonance that should appear in 
real ecosystems, because of geological cause and the environmental noise [HSl- The 
dynamics of the two species is analyzed by fixing the additive noise intensity at 
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Figure 2 . Time evolution of the interaction parameter for differ- 
ent values of the additive noise ap. (a) ap = 0; (b) = 1.78-10"*; 
(c) cr/3 = 1.78 • 10-3; (d) ap = 1.78 • 10"^. The values of the pa- 
rameters are: 7 — 10"^, wo/(27r) ~ 10^'^. 



the value ap — 1.78 • 10"'^, corresponding to a competition regime between the two 
species periodically switched from coexistence to exclusion. The temporal series of 
the two species are obtained for different values of the multiplicative noise intensity 
= (Tx = cFy. The initial values of the two species are x(0) = y(0) = 1. In 
Fig. 3, we report the time series of the two species densities for different values of 
the multiplicative noise. For cr ~ (see Figs. 3a), a regime of coexistence with 
correlated oscillations between the two species is observed. Increasing the intensity 
of the multiplicative noise anti-correlated oscillations appear characterized by a 
larger amplitude with periodical random inversions of populations (see Fig. 3b-3c). 
For higher levels of the multiplicative noise a degradation of the signal and a loss of 
coherence of the temporal series for the species appears (see Fig. 3d). These series 
indicate the presence of stochastic resonance (SR): because of a bistable potential 
modulated by a weak periodic force, the response of the system may be enhanced 
by the presence of the noise and a periodicity appears. We investigate the presence 
of SR by considering {x — y)"^, the squared difference of population densities. Fig. 
4 shows the SNR of this quantity as a function of the multiplicative noise intensity 
cr, for (7p ~ 1.78 • 10"'^. We note that dynamics of {x — y) is mainly affected by 
the multiplicative noise, as we can see from Equations (1), (2). A maximum at 
cr = lO^** is present. The above analysis makes it clear the role of the two noise 
sources: the additive noise determines the conditions for the different dynamical 
regimes of the two species, and the multiplicative noise produces a coherent response 
of the system by a mechanism of symmetry breaking of the dynamical evolution of 
the ecosystem. 
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Figure 3. Time evolution of both populations at different levels of 
the multipHcative noise: (a) a = 0; (h) a = 10~^°; (c) cr = 10^^; 
(d) (7 = 10~^. The values of the parameters are /i = 1, a = 1, 
7 = 10^^, o^o/Stt = 10^'^. The intensity of the additive noise is 
fixed at the value = 1.78 • 10^'^. The initial values of the two 
species are a;(0) = y{0) = 1. 



2.2. Noise-delayed extinction. We now consider the mean extinction time of 
one species as a function of the additive noise intensity ap, by fixing a low value 
of multiplicative noise in such a way that the system is far enough from the SR 
regime ^1 . We are not interested in the coherent behavior of the ecological 
system, but we are focused on the effect of the additive noise on the average ex- 
tinction time of the species. In Fig. 5 the usual initial condition, /3(0) = 0.94, is 
fixed. We note that for ap = the ecosystem is in the coexistence regime; that is, 
the deterministic extinction time of both species is infinite. By introducing noise 
causes exclusion to take place and a finite mean extinction time (MET) appears. By 
varying the intensity of the additive noise in Equation|31we obtain, of course, a vari- 
ation of the average extinction time. The delayed extinction is obtained for noise 
intensities ranging from the intermediate regime (2 in Fig. 5a) to the coexistence 
regime obtained with higher values of a (3 in Fig. 5a). This may mimic the be- 
havior of real ecosystems, where a finite mean extinction time may appear because 
of the presence of a nonvanishing level of noise intensity. For some environmental 
reason the noise intensity can change considerably, as is observed in experimental 
data of populations in a very long time interval |34| . Therefore the dynamical be- 
havior shown in Fig. 5 should explain such physical situations, where the variation 
of the environmental noise produces a delayed extinction of some population. By 
increasing the noise intensity we obtain noise-delayed extinction and the average 
extinction time grows reaching a saturation value, that corresponds to a situation 
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Figure 4. Log-Log plot of SNR as a function of the multiplicative 
noise intensity. The SNR corresponds to the squared difference of 
population densities (x — y)'^. The values of the parameters are the 
same as Figure 13 



in which the potential barrier is absent. We find nonmonotonic behavior of the 
MET as a function of the noise intensity ap, with a minimum value Tmin = 40.47 
at = 2.75 • 10~^, which is of the same order of magnitude of the barrier height h 
(see Fig. 5a). The Kramers time corresponding to this noise intensity is — 41.6, 
which is approximately equal to Tmin- This result is due to the noise-driven dy- 
namics. In fact, for a low value of noise intensities, the average time to overcome 
the potential barrier is very high; that is, Kramers times are long. The ecosystem 
remains in the coexistence regime for a long time and the extinction time is very 
large. For noise intensity of the same order of magnitude of the barrier height, the 
system goes toward the exclusion regime of one of two species, and the average 
extinction time is approximately equal to the Kramers time. We get the minimum 
value of MET. For higher values of noise intensity, the Kramers time becomes very 
small, and the representative point of the (3 parameter moves between the two min- 
ima in a very short time. In this condition the system "sees" the average value of 
the interaction parameter (/? = 0.99), which gives a coexistence regime. In Figs. 
5b, 5c, and 5d we show the time evolution of the ecosystem corresponding to the 
points 1, 2 and 3 of Fig. 5a. We have a coexistence regime in points 1 and 3 and 
an exclusion regime in point 2. 

2.3. Colored noise. In real ecosystems the external random perturbations, be- 
cause of interaction with the environment, are correlated within a finite correlation 
time. When the time scale of random fluctuations is larger than the characteristic 
time scale of the ecosystem the external noise cannot be considered white noise. A 
strongly correlated noise, for example, emerges as the result of a coarse graining 




Figure 5. (a) Mean extinction time of one species as a function of 
the noise intensity a/j. Time evolution of both species for different 
levels of additive noise: (b) = 10~^, (c) CT/j = 2 • 10""^, (d) 
ap = 10^^. The values of the parameters are = 1, a = 45, 
7 = 10"^, ujq/2tt = 10~^. The intensity of the multiplicative noise 
is fixed at the value a = 10~^. The initial values of the two species 
are a;(0) = y(0) = 1. 



over a hidden set of slow variables jJUl- In this section we report the effect of re- 
alistic noise in the dynamics of two competing species, and specifically on the SR 
phenomenon in population dynamics in the presence of exponentially correlated 
noise. The dynamics of our ecosystem is described by Equations (1), (2), and (3). 
For low values of the correlation time Tc, the response of the system coincides with 
that obtained with multiplicative white noise. For higher values of the coherent 
response of the system and the maximum of the signal-to-noise ratio (SNR), which 
are signatures of the SR phenomenon, are shifted towards higher values of the noise 
intensity. These results agree with previous theoretical and experimental investi- 
gations of the SR phenomenon in dynamical systems in the presence of colored 
noise ^3 1351 . However in previous studies the colored noise was additive, while 
here we have two different sources of noise and only one of them is colored. 

Now in Equations (1) and (2), ^i{t) {i = x,y) are colored noises given by the 
archetypal source for colored noise; that is, exponentially correlated processes given 
by Ornstein-Uhlenbeck process [23 

^ = --C. + -mit) {i = x,y) (5) 

at Tc Tc 

and r]i(t) {i = x, y) are Gaussian white noises within the Ito scheme with zero mean 
and correlation fimction {i]i{t)rij{t')) = 2a5{t — t')Sij. The correlation function of 
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the processes of Equation 101 is 

(C«We.(i')) = -e-l*'*'l/^^% (6) 

and gives 25{t — t')5ij in the hmit Tc —t 0. Analogous to the previous case (i.e. 
niultiphcative white noise), the time series for the two populations are obtained 
setting Q!i = a2 = a, = (i^it) = P{t), and £,x{t) = ^yit) = ^{t), where 

the interaction parameter (3{t) is described by Equations (3), (4). The optimum 
coherent time behavior of /3{t) (Fig. 2c), typical of the SR phenomenon, may 
be used to obtain the time series of the two species densities in the presence of 
multiplicative colored noise. Therefore we follow a procedure analogous to that 
applied in the case of multiplicative white noise: we analyze the dynamics of the 
two species by fixing the additive noise intensity at the value ap = 1.78 • 10"'^ (see 
Fig. 2c), and we vary the intensity of the multiplicative colored noise. We obtain 
the time series of the two species for different values both of the multiplicative 
noise intensity a = cFx = <^y and the correlation time Tc |14I I38| . In particular 
we investigate the system for (a) Tc < To and (b) Tc > To, with To the period 
of the deterministic driving force. In the weak correlated noise regime (a), no 
relevant modifications occur in the temporal series of the two species densities in 
comparison with the case of multiplicative white noise. The time evolution of 
the two species shows an anticorrelated behavior with quasiperiodical oscillations 
with a random inversion of the population that predominates over the other one, 
as in the white noise case. For Tc ~ To some modifications occur. In particular 
for Tc = 2 • 10'^ the time series of the two species densities show anticorrelated 
behavior with quasiperiodical oscillations up to ct = 10~^ (see Fig. 6d); that is, a 
delay in the coherent output of our ecosystem. This delay will manifest itself in 
the behavior of the signal-to-noise ratio (SNR) as a function of the multiplicative 
noise intensity. In the strong correlated noise regime (b), a relevant delay of the 
coherent time behavior of the two species is observed. The maximum SNR is 
shifted toward higher values of the multiplicative noise intensity. This shift in a 
Log-Log scale grows faster than a linear function of the correlation time Tc- The 
coexistence regime and the correlated oscillations of both populations persist for 
a wider range of multiplicative noise intensities. The anticorrelated behavior with 
quasiperiodical oscillations appears with very high noise intensity as the correlation 
time value of the multiplicative noise is strong enough. The loss of coherence in 
the time behaviors of the two species is observed at very high intensities of the 
multiplicative noise. Because of the high values of the multiplicative noise, one 
population extinguishes and the other one survives at a constant density after a 
transient dynamics. This dynamical behavior is typical of an ecosystem in the 
presence of an absorbing barrier According the case of multiplicative noise, 
to underline the presence of SR, we analyze the squared difference of population 
densities {x — y)^ for different values of Tc. In Fig. 7 the SNRs of this quantity 
are shown for = 0, 2 • 10-^,2 • 10-\20,2 • 10^,2 • 10^ 2 • 10^ as a function of 
the multiplicative noise intensity a, by fixing the additive noise intensity |38j at 
ap = 1.78 • 10~^. In each graph of this figure a maximum appears, whose position 
depends on the values of Tc', that is, the most coherent response of the system 
is connected with both the intensity and the correlation time of the multiplicative 
noise. We see clearly the two dynamical regimes: (a) weak correlated noise (the first 
four values of Tc), (b) strong correlated noise (the last three values of Tc). In this 
second regime the maximum SNR is shifted toward higher values of multiplicative 
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Figure 6. Time evolution of both populations at different levels of 
the multiplicative noise for Tc — 2 ■ 10'^: (a) (7 = 0; (b) tr = 10~^^; 
(c) a = 10~^; (d) a ~ 10~^. The values of the parameters are 
= 1, a = 1, 7 = lO^S wo/(27r) = 10"^. The intensity of the 
additive noise is fixed at the value ap = 1.78 • 10^'^. The initial 
values are: for the two species a;(0) = y(0) = 1, for the additive 
(white) noise /3(0) = 0.94, for the multiplicative (colored) noise 

Ci(o) = C2(o) = o. 



noise intensity as in previous theoretical and experimental studies |10[ 1111 I36j . 

However some differences occur. Previous studies on the effect of colored noise on 
the SR phenomenon showed that by increasing Tc the peak of the SNR shifts toward 
higher values of the noise amplitude and the maximum decreases with a broadening 
of the entire curve. The shift of the SR peak to larger noise intensities is due to the 
fact that colored noise suppresses exponentially the hopping rate with increasing 
noise color. In our model the colored noise is introduced in the multiplicative noise 
and not in the additive one as in usual bistable dynamical systems. The SR in the 
dynamics of the interaction parameter (3 induces SR phenomenon in the dynamics of 
two competing populations |14II38| . Our hopping rate in the first SR is not affected 
by the "color" of the multiplicative noise. However, this noise is responsible for 
the coherent response of the ecosystem, and therefore the presence of color in the 
multiplicative noise causes the SNR peak to shift. 



2.4. Spatially extended systems. To study the spatial effects due to the pres- 
ence of noise sources we consider a discrete time evolution model, which is the 
discrete version of the Lotka-Volterra equations with diffusive terms, namely a cou- 
pled map lattice (CML) ^ 
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Figure 7. Log-Log plot of SNR as a function of noise intensity. 
The SNR is obtained for six different values of the correlation time: 
Te = 2 • 10-2, = 2 • 10-1; Tc = 20; r, = 2 • 10^, r, = 2 • 10^; 
Tc = 2 • 10^. Moreover the signal-noise ratio for white Gaussian 
noise is reported. The SNR corresponds to the squared difference 
of population densities (x — y)^. The values of the parameters and 
the initial condirtions are the same of those used to obtain the 
temporal series. 



7 

y.";' = nyl^{l-yl^-l3-xl^) + ^yyl^Y::^+DY,{y-~y'l^). (8) 

7 

In Equations ^ and (O x"^ and y"^- denote respectively the densities of species x 
and species y in the site at the time step n, ^ is proportional to the growth 
rate, D is the diffusion constant, indicates the sum over the four nearest neigh- 
bors. The random terms are white noise sources, modeled by independent Gaussian 
variables denoted by Xf^, Y"j with zero mean and variance unit. Here ax, <Jy are 
the intensities of the multiplicative noise that models the interaction between the 
species and the environment. The interaction parameter /3" of Equations Q) and 
(jSJ is a stochastic process which corresponds to the value of continuous of 
Equation taken at the step n and u)q/2tt = 10-^. 

We consider the time evolution of the spatial distribution of the ecosystem, 
described by Equations ((JJ) and (jHl , in the SR dynamical regime obtained for ap = 
2.65 • IQ-^. We fix the additive noise at this value and vary the intensities of 
multiplicative noise. 



12 



B. SPAGNOLO, D. VALENTI, A. FIASCONARO 



We obtained spatio-temporal patterns of the two species for different values of the 
multiplicative noise intensity a = ax = (Ty, namely a = 10~^^, 10~^, 10""*, 10^^ with 
fi = 2, D = 0.05, 7 = 1.5 • 10-1, wo/(27r) = lO^^, /3(0) = 0.94 and = yl^ = 0.5 
at all sites )40| . For very low noise intensity an average correlation on the 

considered lattice {N = 100 x 100) between the species is observed. For higher noise 
intensities an anticorrelation between the two species is observed: the two species 
tend to occupy different positions. The anticorrelation is more evident for a = 
lO"'^ (see Fig. 8a). Increasing the multiphcative noise reduces the anticorrelation 
strongly (see Fig. 8b). 
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Figure 8. Spatial distributions at different times for (a) a = 10^^ 
and (b) a = 10"^. The value of the additive noise is fixed at = 
2.65 • 10^'^. The values of the parameters are: /i = 2, D = 0.05, 
7 = 1.5 • 10-\ wo/(27r) = 10-2, iV = 100 x 100. The initial values 
are x^^ = = 0.5 for all sites and /3(0) = 0.94. 



To evaluate the spatial correlation between the two species for the noise intensi- 
ties considered, we calculate, at the time step n, the correlation coefficient < c" > 
defined on the lattice as 001 



< c" (9) 



with 



coti„ 



E..,-«j-x")(yi:,--r) 

N 



(10) 



where i", s", y", are the mean value and the root mean square respectively of 
species 1 and species 2, obtained over the whole spatial grid at the time step n, 
cow" is the corresponding covariance and N = 100 x 100 the number of sites which 
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compose the grid. The behavior of the correlation coefficient < c" > as a function 
of the time for different levels of the multiplicative noise is reported in Fig. 9 |4(J| . 
We observe a nonmonotonic behavior of < c" > as a function of the multiplicative 
noise intensity. For low noise intensities cr — 10^^^, < c" > shows weak oscillation 
around 1, that is strong correlation between the two species. For higher levels 
of the noise cr = 10^^", < c" > is affected by fluctuations and its values vary 
strongly as a function of the time. A further increase of the multiplicative noise, 
(i.e., cr = 10~^ and cr = 10"''), determines an oscillation of < c" > around a 
negative value (i.e., anticorrelation between the two species), with the frequency of 
the periodical forcing. For higher intensities of the noise cr = 10"^, the value of 
the correlation coefficient < c" > increases and it vanishes for a = 10+'^. To show 
clearly the nonmonotonic behavior of < c" > , we calculate the time average of the 
correlation coefficient < c" >t and we report it, as a function of the multiplicative 
noise intensity, in Fig. 10. A clear minimum appears, which corresponds to the 
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Figure 9. Correlation coefficient < c" > as a function of the time. 
For low levels of the multiplicative noise (cr = 10^^^) the species 
are strongly correlated and < c" > is approximately constant. 
By increasing the intensity of the multiplicative noise (cr ~ 10^^") 
< c" > shows big fluctuations. A further increase of the noise (cr ~ 
10~^, a = 10~^) causes strong anticorrelation between the two 
species with < c" > oscillating at the frequency of the periodical 
forcing. For very high levels of noise, the anticorrelation is reduced 
(cr = 10~^) and finally it disappears (cr = 10+'^); that is, the species 
are totally uncorrelated. 

anticorrelated oscillations shown in the time evolution of two competing species in 
each point of our spatial grid. We note therefore the different role of the two noise 
sources in the ecosystem dynamics. The additive noise determines the conditions of 
the dynamical regime, the multiplicative noise produces a coherent response of the 
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system ^1 OHl , which is responsible for the appearance of anticorrelation behavior 
in the spatial patterns of the species. 




3. Three interacting species. In this section we report the spatio-temporal dy- 
namics of three interacting species, two preys and one predator, in the presence of 
multiplicative white noise and a periodical driving force. We use the same coupled 
map lattice model of the previous section 

<r = - - /3"y", - 1^1,) + + D - (11) 

s 

where x"^ , y"^, and zf^ are, respectively, the densities of preys x, y and of the 
predator z in the site at the time steps n. Here 7 and 7^ are the interac- 
tion parameters between preys and predator and D is the diffusion coefficient. In 
previous equations X , Y and Z are the white Gaussian noise variables with 

{Xit)) = {Yit)) = {Z{t)) = 0, (14) 

{X(t)X{t + r)) = {Y{t)Y{t + r)) = {Z{t)Z{t + r)) = <5(r), (15) 

<^x = cTy = = g is the noise intensity, and /i and /i^ are scale factors. indicates 
the sum over the four nearest neighbors in the map lattice. The boundary conditions 
have been established in such a way that no interaction is present out of lattice. 
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The interaction parameter (3 between the two preys is a periodical function whose 
value, after n time steps, is given by 

f3{t) = l + e + acos{ujot), (16) 
with e — —0.01, a — 0.1, and — {10^/211) — IQ^^. The interaction parameter 




prej" X p roy y prcd s.toc z 



Figure 11. Spatial patterns induced by the noise for three in- 
teracting species (two preys and one predator) with homogeneous 
initial distributions. The parameter set is: e = —0.01, fi = 2, 
fi, = 1,1^ = 1, (3, = 0.01, z/o = (c^o/27r) = IQ-^, a = 0.1, a, - 
ay=a. = 10"^ D = 0.01, 7 = 3- 10"^, 7^ = 2.05 • 10^. The initial 
values of the uniform spatial distribution are x^™* = y*™* = 0.25 
and z^f = 0.10 for all sites 

(3(t) oscillates around the critical value /3c = 1 in such a way that the dynamical 
regime of Lotka-Volterra model for two competing species changes from coexistence 
of the two preys (/3 < 1) to exclusion of one of them {(3 > 1). We consider two 
different initial conditions: (i) a homogeneous initial distribution, and (ii) a peaked 
initial distribution. In the first case we find exactly anticorrelated spatial patterns 
of the two preys, while the spatial patterns of the predator show correlations with 
both the spatial distributions of the preys (see Fig. 11). The preys tend to occupy 
different positions as in the case of two competing species. In the second case, we 
use delta-like initial distributions for the two preys and a homogeneous distribution 
for the predator. After 800 steps we find strongly correlated spatial patterns of 
the preys that almost overlap each other. The maximum of spatial distribution of 
the predator is just at the boundary of the spatial concentrations of the preys, so 
that the predator surrounds the preys (see Fig. 12). The preys now tend to overlap 
spatially, as occurs in real ecosystems when preys tend to defend themselves against 
the predator attacks [Oj . 

The quantitative calculations of the site correlation coefhcient between a couple 
of species in the lattice have been done using the following formula 

^ EZ(^",~^"){K,~kn ^ ^^^^ 

where TV is the number of sites in the grid, the symbols w", fc" represent one of the 
three species x, y, z, and ?z}". A:" represent the mean values of the concentration of the 
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prey X prey y pire dai or z 



Figure 12. Spatial patterns induced by the noise for three inter- 
acting species (two preys and one predator) with delta-hke initial 
distributions of the preys and a homogeneous distribution of the 
predator: (a) initial conditions, (b) spatial patterns after 800 time 
steps. Here we set e = —0.05, D — 0.1, (Tx ^ (Ty = (Tz = lO^'' and 
the other parameters are the same as in Figure [TTl 

species in all the lattice at the step n. The two-dimensional spatial grid considered 
is composed by = 100 x 100 sites in {x,y) plane. The calculations have been 
done for various noise intensities and at different steps of the iteration process. 
To quantify our analysis, we consider only the maximum patterns, defined as the 
ensemble of adjoining sites in the lattice for which the density of the species belongs 
to the interval [3/4 max, max], where max is the absolute maximum of density in 
the specific grid |41| . For each spatial distribution, in a temporal step and for a 
given noise intensity value, the following quantities have been evaluated referring 
to the maximum pattern (MP): mean area of the various MPs found in the lattice 
and spatial correlation r between two preys, and between preys and predator. The 
parameters used in our simulations are as follows: a = 0.2, ujq = ttIO^'', e = —0.1 
= 2, = 1, 7 = 0.03, fiz = 0.02, 7^ = 205, and D = 0.1. The noise intensity 
<^x — o'y — o'z varies between 10""'^^ and 10~^. 

3.1. Deterministic analysis. In the absence of noise and setting constant the 
value of the interaction parameter f3 we obtain: (i) for e < (/3 < 1) a coexistence 
regime of the two preys characterized in the lattice by a strong correlation between 
them with the predator lightly anticorrelated with the two preys; (ii) for e > 
(/3 > 1) wide exclusion zones in the lattice (see Fig. 13), characterized by a strong 
anticorrelation between preys. 

By considering the periodic variation of the interaction parameter /?(t), we ob- 
tain for e = 0, after a transient anticorrelated behavior between preys, a coexistence 
regime with strong correlation between preys that evolves toward an homogeneous 
spatial distribution of all three species. For e > 0, we find an oscillating behavior 
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Figure 13. Site correlation coefficient r in noiseless dynamics 
as a function of time for different values of the parameter e : 
— 0.1, 0., +0.1. Here rj = 0.2. The parameter set is: P = 1.1,(7 = 
0,L» = 0.1, yu = 2, 1/ = l,a = 0.03, fi^ = 0.02,7 = 205. The initial 
conditions arc random with a Gaussian distribution, with mean 
values x(0) = y{0) = z{0) — 0.25 and variance ctq = 0.1. Here 
^12, i~i3,r23 and ri23 are respectively the site correlations between: 
(i) preys, (ii) prey 1 and predator, (iii) prey 2 and predator, and 
(iv) predator and both preys. 



of the site correlation coefficient from coexistence regime between preys (corre- 
sponding to strong correlation) to an exclusion regime, corresponding to strong 
anticorrelation. This last behavior is prevalent. The oscillating frequency coincides 
with that of the /3-parameter. When e < 0, the two preys, after an initial transient, 
remain strongly correlated for the entire time, despite the fact that the parameter 
(3(t) takes values greater than 1 during the periodical evolution. This situation 
corresponds to a coexistence regime between preys. In Fig. 13, we report the be- 
havior of the site correlation coefficient r as a function of time for three values of 
the parameter e = —0.1,0,0.1 [411 . 

3.2. Spatial patterns induced by noise. To analyze the effect of the noise, 
we focus on the interesting dynamical regime characterized, in absence of noise, 
by coexistence between preys in all the period of /? that is, with e < 0. The 
noise triggers the oscillating behavior of the site correlation coefficient r giving rise 
to periodical alternation of coexistence and exclusion regime. Even a very small 
amount of noise is able to destroy the coexistence regime periodically. Noise is 
also responsible for a nonmonotonic behavior of the area of spatial patterns, which 
repeats periodically in time. In Fig. 14, we report a nonmonotonic behavior of the 
area of the maximum pattern as a function of noise intensity. A maximum of the 
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Figure 14. Semi-log plot of the mean area of the maximum 
patterns for all species as a function of noise intensity, at itera- 
tion step 1400. Here circles and triangles are related to preys, 
squares to predator and e = — 0.1,77 = 0.2. The values of the 
other parameters are the same used for Figure 13. The initial 
spatial distribution is homogeneous and equal for all species, i.e. 
X™* = y™* = z™* = 0.25 for aU sites 



area of maximum patterns is visible for the preys aX q ~ 10^^ and for the predator 
at q = 10^^. The same behavior is present in the following time steps within the 
first period of the interaction parameter: 600, 800, 1200, and 1400. But at time 
steps 600 and 800 the preys are highly uncorrelated with site correlation coefficient 
ri2 = — 1, while at time steps 1200 and 1400, the preys are highly correlated with 
ri2 = 1 . The formation of spatial patterns appears only when the preys are highly 
correlated, while large patches with clusterization of preys appear when they are 
uncorrelated. This means that the coexistence regime between preys corresponds 
to the appearance of spatial patterns, while the exclusion regime corresponds to 
clusterization of preys. The noise- induced pattern formation relative to the iteration 
1400 is visible in Fig. 15, where we report five patterns of one prey and the predator 
for the following values of noise intensity: q = 10~^^, 10~^, 10~*, 10"'', 10~^. The 
initial spatial distribution is homogeneous and equal for all species; that is, a;^™* = 
y'lmt _ ^mit _ Q 25 foi- all sites A spatial structure emerges with increasing 

noise intensity. This spatial pattern disappears for sufficiently large noise intensity 
(see Fig. 15e). As a final investigation, we analyze the behavior of the area of 
the patterns as a function of time. We observe a nonmonotonic behavior of the 
area of MPs as a function of time for all values of the noise intensity investigated. 
Particularly for noise intensity values greater than q = 10^^, this nonmonotonic 
behavior becomes periodic in time with the same period of P{t), as shown in Fig. 
16 for q = 10~^. We note that this nonmonotonic behavior doesn't mean that a 
spatial pattern appears, like that of Fig. 15b, but that a big clusterization of preys 
density may occur. The maximum at g = 10~^ of Fig. 16a in fact corresponds to 
large patches of preys in the lattice investigated 01]. The various quantities, such 
as pattern area and correlation coefficient, have been averaged over 200 realizations, 
obtaining the mean values shown in the Figs. 14 and 16. The effects induced by the 
multiplicative noise can be summarized as follows: (i) to break the symmetry of the 
coexistence regime between the preys, producing an alternation with the exclusion 
regime; (ii) to trigger the oscillating behavior of the site correlation coefficient; and 
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Figure 15. Spatial pattern formation for preys and predator, at 
time iteration 1400 and for the following values of the noise inten- 
sity: q = 10"", 10-3, 10-^ 10"'', 10-2. The values of the other 
parameters and the homogeneous initial distribution are the same 
used in Figure 14. The parameters ri2, ri^, r23, ri23 have the same 
meaning of Figure 13. 



(iii) to produce a nonmonotonic behavior of the pattern area as a function of the 
noise intensity with an appearance of spatial patterns. 

4. N interacting species. In the last part of this short review we report the main 
results obtained by analyzing an ecosystem composed by N interacting species in 
a noisy environment in the presence of an absorbing barrier; that is, extinction of 
the species PdlH]. 
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Figure 16. Mean area of Maximum pattern of the three species 
and relative sites correlations between preys and between preys 
and predator as a function of time and q = 10~*. (a): black 
and white circles are related to preys, triangles to predator; (b) 
site correlation coefficient ri2 (black circles), ri3 (white circles), 
and ri23 (triangles). The values of the other parameters and the 
homogeneous initial distribution are the same used in Figure 13. 
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We consider an N-species generalization of the usual Lotka-Volterra system, and 
the Ito stochastic differential equation describing the dynamical evolution of the 
ecosystem is 



ni{t)dt + ^/eni{t)dwi, i=l,...,N (18) 

where ni(i) > is the number of elements of the i*'' species. In Equation (|18() . 7 is 
the growth parameter, the interaction matrix Jy modelizes the interaction between 
different species (i ^ j), and Wi is the Wiener process whose increment dwi satisfies 
the usual statistical properties 

< dwi{t) > = 0; < dwi(t)dwj(t') > ^ 5.,j5{t - t')dt. (19) 



dn^{t) 



n,(t) +^ Jynj(i) 



We consider all species equivalent so that the characteristic parameters of the 
ecosystem are independent of the species. The random interaction with the environ- 
ment (climate, disease, etc..) is taken into account by introducing a multiplicative 
noise in the Equation H18|) . The solution of the dynamical Equation (|18|) is given 

by 



ni{0)exp 


5t^ 






1 + 771^(0) /jj* dt'exp 


"<5t' + Viu,,(t')+/o' rfi"E,^.^.. 





(20) 



We consider two different types of interaction between the species: (a) a mean 
field approximation with a symbiotic interaction between the species; (b) a ran- 
dom interaction between the species with different types of mutual interactions: 
competitive, symbiotic, and prey-predator relationship. 
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4.1. Mean field approximation. We consider a mean field symbiotic interaction 
between the species. As a consequence, the growth parameter is proportional to 
the average species concentration, 



In the limit of a large number of interacting species the stochastic evolution of the 
system is given by the integral equation 



M{t) = 1 Vln fl + n,(0) /* dt'e'^A^(*')+7*'+V^».(t') 



(22) 



where M{t) is the time integral of the site population concentration average. We 
introduce an approximation of this Equation l|22() which greatly simplifies the noise 
affected evolution of the system and allows us to obtain analytical results for the 
population dynamics. We note that in this approximation the noise influence is 
taken into account in a nonperturbative way, and that the statistical properties of 
the time average process M(t) are determined asymptotically from the statistical 
properties of the process Wmax{t) = supo<t/<t w(t'), where w is the Wiener process. 
Starting from the following approximated integral equation for M(t) 



M{t) 



l + n,{0)e 



•/ewr, 



(23) 



it is possible to analyze the role of the noise on the stability-instability transition 
in three different regimes of the nonlinear relaxation of the system: (i) toward the 
equilibrium population (7 > 0); (ii) toward the absorbing barrier (7 < 0); (iii) at 
the critical point (7 = 0). Specifically at the critical point we obtain for the time 
average process M{t)/t as a dominant asymptotic behavior in the stability region 
(namely when J < 1) 



M{t) _ _ 
and in the instability region (namely when J > 1) 




(24) 



Af(t) 



= (ln(n,(0))> 



277 e 



t 



Vi 



Vt 



(25) 



4.2. Random interaction. The interaction between the species is assumed to be 
random and is described by a random interaction matrix J^, whose elements are 
independently distributed according to a Gaussian distribution 



1 



where J is the interaction strength and 



exp 



N 



(26) 



< Ji, 0, 



<C JijJji 



>= 0. 



(27) 



With this asymmetric interaction matrix, our ecosystem contains 50% prey-predator 
interactions (namely Jy < and Jji > 0), 25% competitive interactions (Jy < 
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and Jji < 0), and 25% symbiotic interactions (J^ > and Jji > 0). The initial 
values of the populations ni(0) have also Gaussian distribution 



P{n) 



zcxp 



(n- < n >)^ 



2ai 



0.01, and < n >= 1. 



(28) 



The strength of interaction between the species J determines two different dynam- 
ical behaviors of the ecosystem. Above a critical value Jc, the system is unstable; 
this means that at least one of the populations diverges. Below the critical interac- 
tion strength, the system is stable and asymptotically reaches an equilibrium state. 
For our ecosystem this critical value is approximately J = 1.1. The equilibrium 
values of the populations depend both on their initial values and on the interaction 
matrix. If we consider a quenched random interaction matrix, the ecosystem has a 
great number of equilibrium configurations, each one with its attraction basin. For 
vanishing noise (e = 0), the steady state solutions of Equation H18|l are obtained by 
the fixed-point equation 



(7 - ni + hi)nt ^ 



(29) 



where 



/li = ^ Jijnj{t) 



(30) 



is the local field. For a large number of interacting species, we can assume that the 
local field hi is Gaussian with zero mean and variance a\. =< hj >= J"^ < nf > 



The solutions of Equation (|29|) are 



^exp 



2< 



(31) 



and 



rii = 0, i. e. extinction 



rii = + hi)0{"/ + hi), Hi > 0, 



(32) 
(33) 



where Q is the Heaviside unit step function. From this equation and applying the 
self-consistent condition, we can calculate the steady state average population and 
its variance. Specifically, we have 



<n,> = {{j + h,)e{-f + h,)) ^ 




,(34) 



and 
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For an interaction strength J ~ I and an intrinsic growth parameter 7 = 1 we 
obtain: < >= 1.4387, < nf >= 4.514, and a^. = 2.44. These values agree 
with that obtained from numerical simulation of Equation l|18ll . The choice of this 
particular value for the interaction strength, based on a preliminary investigation on 
the stability-instability transition of the ecosystem, ensures us that the ecosystem 
is stable. The stationary probability distribution of the populations is the sum of 
a delta function and a truncated Gaussian 



P{ni) = ne,5{ni) + 6(71^)- 



exp 



(?I.-Tlio)^ 
2J2(j2 



(36) 



The stationary probability distribution of the population densities has been ob- 
tained, without the extinct species, in comparison with the computer simulations 
for systems with N = 1000 species and for an interaction strength J = 1, and for 

7 = 1 IHl- 

Now we focus on the statistical properties of the time integral of the i*'' popula- 
tion Ni{t) 



dt'n,{t'), 



(37) 



in the asymptotic regime. From Equation H20|l we have 



iV,(i) =ln 



1 + nj(0) / dt'exp 



7i' + \/im,(t') + ^ Jy^j(t') 



(38) 



In Equation (??) the term J^j Jij^j gives the influence of other species on the 
differential growth rate of the time integral of the i*'* population and represents a 
local field acting on the i*'' population ^ |S1 031 



(39) 



We use the same approximation of the Equation (|23|l and, after differentiating, we 
get the asymptotic solution of Equation (??) 



NS) - In 



(40) 



where Wmax.it) = supo<t' <tw{t') and r]m.ax,it) = sMpo<t'<t?7(i')- Equation (gUJ) is 
valid for 7 > 0; that is, when the system relaxes toward an equilibrium population 
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and at the critical point. Evaluating Equation Ij4()|l for 7 > 0, after making the 
ensemble average, wc obtain for the time average of the i*^ population Ni 



and 



Ny_,Vet+{"f + Nr, + (ln 



^{o)]) 


, 7 


ni{o) 


)] 


. 7 . 





, 7>0, 



(41) 
(42) 



where and Njj are variables with a semi-Gaussian distribution ^ and Njj must 
be determined self-consistently from the Equation H39|) . 

These asymptotic behaviors are consistent with those obtained using a mean 
field approximation. We obtain in fact the typical long time tail behavior {t^^^^) 
dependence, which characterizes nonlinear relaxation regimes when 7 > 0. Further, 
the numerical results confirm these analytical asymptotic behaviors of Ni |27| . 

When the system relaxes toward the absorbing barrier (7 < 0) we get from 
Equation (??) in the long-time regime 



if.) = f 



ln(n,;(0)) 4-ln 



dt'e 



(43) 



In this case the time average of the i*'' population {Ni) is a functional of the local 
field and the Wiener process, and it depends on the history of these two stochastic 
processes. We have also analyzed the dynamics of the ecosystem when one species 
is absent. Specifically, we considered the cavity field, which is the field acting on the 
i*'' population when this population is absent |43j . The probability distributions 
for both local and cavity field have been obtained by simulations for a time < = 100 
(expressed in arbitrary units) in absence of external noise, and for two species 
(namely species 1 and 33) 0. We found that the probability distributions of the 
cavity fields differ substantially from those of local fields for the same species unlike 
the spin glasses dynamics, where the two fields coincide. The same quantities have 
also been calculated in the presence of the external noise |H|. The effect of the 
external noise is that the two fields overlap in such a way that for some particular 
species they coincide. This interesting phenomenon, which is reminiscent of the 
phase-transition phenomenon, was found for some populations. The main reasons 
for this peculiar behavior are (i) all the populations are positive; (ii) the particular 
structure of the attraction basins of our ecosystem; and (iii) the initial conditions, 
which differ for the value of one population, belong to different attraction basins. 
Some populations, in the absence of external noise, have a dynamical behavior such 
that after a long time they significantly influence the dynamics of other species. In 
the presence of noise, all the populations seem to be equivalent from the dynamical 
point of view. We also found that for strong noise intensity (namely e = 1) all 
species extinguish on a long-time scale [t k, 10^ a. u.). Whether extinction occurs 
for any value of noise intensity is still an open question, because of time-consuming 
numerical calculations. 



5. Conclusions. We briefly reviewed the noise-induced phenomena in population 
dynamics of three different ecosystems: (i) two competing species, (ii) three inter- 
acting species, and (iii) N-interacting species. In the case of two competing species, 
we considered two sources of white noise: a multiplicative noise and an additive 
noise, that produces a random interaction parameter between the species. The 
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noise induces a coherent time behavior of two species, giving rise to temporal os- 
cillations and enhancement of the response of the system through the stochastic 
resonance phenomenon. Specifically the additive noise controls the switching be- 
tween the coexistence and the exclusion dynamical regimes, the multiplicative noise 
is responsible for coherent oscillations of the two species. The SR in the dynam- 
ics of interaction parameter (3 induces SR phenomenon in two competing species. 
These time behaviors are absent in the deterministic dynamics. The noise is also 
responsible for a delayed extinction that gives rise to a nonmonotonic behavior of 
the average extinction time as a function of the additive noise intensity. We evalu- 
ated the role of colored noise and its effects on the time behavior of the two species. 
We found that the multiplicative noise is responsible for periodical oscillations of 
the two species densities, whose amplitude and coherence depend on the value of 
the correlation time Tc- For Tc — > our results are consistent with that obtained 
for the case of white noise. Moreover the coherent time behavior of our ecosystem 
and the SR phenomenon are shifted towards higher noise intensities, in agreement 
with previous theoretical and experimental investigations. We note that our model 
is useful to describe physical situations in which the amplitude of the periodical 
driving force, because of the temperature variations, is weak and therefore unable 
to produce considerable variations of the dynamical regime of the ecosystem. The 
synergetic cooperation between the nonlinearity of the system and the random and 
periodical environmental driving forces produces, therefore, a coherent time behav- 
ior of the ecosystem investigated. We find that these noise-induced effects should be 
useful to explain the time evolution of species whose dynamics are strongly affected 
by the noisy environment |S1 El I42L 151] . We also analyzed the role of the noise in 
spatio-temporal behaviors by using a discrete version of Lotka-Volterra equations 
with diffusive terms. We found that the noise induces spatio-temporal behaviors 
that are absent in the deterministic dynamics; that is, pattern formation with the 
same periodicity of the deterministic force. Moreover appearance of temporal os- 
cillation is observed in the correlation coefficient between the two species and a 
nonmonotonic behavior of the time-average correlation coefficient as a function of 
the multiplicative noise. 

We also analyzed the role of the noise on the spatio-temporal behaviors of an 
ecosystem composed by three interacting species. We found that the formation 
of dynamical spatial patterns occurs with correlations which are strongly depen- 
dent on the initial conditions. Moreover, we obtain nonmonotonic behavior of the 
mean area of the maximum patterns as a function of noise intensity. We find the 
same behavior for the area of the patterns as a function of evolution time. The 
noise changes the dynamical regime of the species, breaking the symmetry of the 
coexistence regime. In addition, the noise produces spatial patterns and temporal 
oscillations of the site correlation coefficient defined on the lattice. Our model for 
spatially extended systems composed by two and three species could be useful to ex- 
plain spatio-temporal behaviors of populations whose dynamics is strongly affected 
by the noise and by the environmental physical variables; that is, interpreting the 
experimental data of population dynamics strongly affected by the noise O 021 ■ 
Finally we analyzed the nonlinear relaxation of an ecosystem composed by N inter- 
acting species. By using an approximation of the integral equation, which gives the 
stochastic evolution of the system, we obtained analytical results that reproduce 
very well almost all the transient dynamics. We investigated the role of the noise 
on the stability-instability transition and on the transient dynamics. For random 
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interaction, we obtained asymptotic behavior for three different nonhnear relax- 
ation regimes. We obtain the stationary probabiUty distribution of the population, 
which is the sum of two contributions: (i) a delta function around n ~ for the 
extinct species and (ii) a truncated Gaussian for the alive species. When we switch 
on the external noise, an interesting phenomenon is observed: the local and the 
cavity fields, whose probability distributions are different in the absence of noise, 
coincide for some populations. This phenomenon can be ascribed to the peculiarity 
of the attraction basins of our ecosystem. This model could be useful to describe 
plankton dynamics. 
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